Well-posedness of some initial-boundary-value problems for 
dynamo-generated poloidal magnetic fields 

[corrected version] 

Ralf Kaiser^ and Hannes Uecker'^ 

^ Fakultat fiir Mathematik, Physik und Informatik, Universitat Bayreuth, D-95440 Bayreuth, 

Germany 

^ Institut fiir Mathematik, CvO Universitat Oldenburg, D-26111 Oldenburg, Germany 
ralf.kaiser@uni-bayreuth.de, hannes. ueckerQuni-oldenburg.de 
December, 2012 

Abstract 

Given a bounded domain G C Mf^, d > 3, we study smooth solutions of a linear 
parabolic equation with non-constant coefficients in G, which at the boundary have 
to C-'^-match with some harmonic function in Mf^ \ G vanishing at spatial infinity. 

This problem arises in the framework of magnetohydrodynamics if certain 
dynamo-generated magnetic fields are considered: For example, in the case of 
axisymmetry or for non-radial flow fields, the poloidal scalar of the magnetic field 
solves the above problem. 

We first investigate the Poisson problem in G with the above described bound- 
ary condition as well as the associated eigenvalue problem and prove the existence 
of smooth solutions. As a by-product we obtain the completeness of the well- 
known poloidal "free decay modes" in M'^ if G is a ball. Smooth solutions of the 
evolution problem are then obtained by Galerkin approximation based on these 
eigenfunctions . 
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1 Introduction 

We are concerned in this paper with the following initial-boundary-value problem: 



dfU — aAu = b ■ Vu + cu in G x M+, (1.1a) 

Au = in G X M+, (1.1b) 

u and Vti continuous in M'^ x M+, (1-lc) 

u{x,t) — )• Uoo{t) for |x| — )■ oo, t G (1-ld) 

ti(-, 0) = no on Gx{t = 0}. (l.le) 



Here, G C M"^, d > 3 is a bounded domain with (sufficiently) smooth boundary dG, and 
G := M'^ \ G. The scalar-valued coefficients a and c, and the vector-valued coefficient b 
are sufficiently smooth functions of x G G and t € M+; a is, moreover, bounded from 
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below by ao > 0. The asymptotic behaviour of solutions at spatial infinity is described 
by the (given) function : — )• M, and the initial-value uq is prescribed on G only. 

Problem (jl.ip arises in the context of magnetohydrodynamic dynamo theory: the 
generation of a magnetic field B by motion of a liquid conductor (of magnetic diffusivity 
r/ > 0) according to some prescribed flow field v is described by the induction equation 

(cf. m) 

dtB = V X {v X B) -V X {r]V X B), V ■ B = 0. (1.2) 

Equation (|1.2p i constitutes a system of parabolic equations for the magnetic field com- 
ponents. In general, the flow field couples these components in a nontrivial way which 
makes the question for "dynamo solutions", i.e. solutions which do not decay in time, 
difficult to answer. Only in special situations does a field component or a related scalar 
quantity decouple, and a general decay result, a so-called antidynamo theorem, may 
be obtained. For instance, if the conductor fills a ball Br C with radius R, if the 
conductivity is radially symmetric, and if the flow field has no radial component the 
quantity P := B ■ x satisfies the scalar problem (cf. appendix A or |llj): 



dtP - r]AP = -V • (vP) 


in Br X 


M+, 


(1.3a) 


AP = 


in Br X 


M+, 


(1.3b) 


P and VP continuous 


in X I 




(1.3c) 


P{x,t) = 0{\x\-^) 


for X — > 


■ oo, t G 


(1.3d) 


P(-,0) = Po, {Po) = 


on Br X 


{t = 0}. 


(1.3e) 



Note that conditions (jl.lb ) and (|l.m ) with Uoo = imply the spatial decay condition 
P{x,t) = 0{\x\-^) (see Appendix C). The stronger condition ()1.3d ) is a consequence 
of the additional zero-spherical-mean condition 

{Po):=^ Pods = 0, re(0,P), 

on the initial value. This condition is preserved by eqs. ()1.3b .b) and holds, consequently, 
for P on M'^ X Equation ()1.3b ) describes a vacuum field outside the conductor and 
condition (|1.3b ) guarantees a continuous magnetic field throughout space. 

Another instance is the axisymmetric dynamo problem in ordinary space M'^. Again, 
a scalar quantity describing the poloidal part of the magnetic field decouples and, if 
reformulated in M^, is precisely a solution of problem (|1.1|) (see appendix B). Note 
that the conductor is here assumed to be axisymmetric but need not be a ball. It 
is this application which motivates the investigation of problem (jl.ip in more than 3 
dimensions and in domains more general than balls. 

The focus of dynamo theory is less on existence theorems than on decay results for 
the magnetic field under certain restrictions on the magnetic field and/or the flow field, 
thus excluding dynamo action under these restrictions. However, proving decay results 
requires sometimes the solution of an auxiliary problem. For instance, in proving a 
"non-radial velocity theorem" for solutions of problem (jl.3p one needs positive solutions 
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of an auxiliary problem of type (jl.ip : and it is this application which requires a non- 
zero asymptotic condition like (jl.ld ) (cf. [llj). Similarly in the axisymmetric problem, 
Backus makes use of solutions of an auxiliary problem (cf . [3] ) . He made the existence of 
such solutions plausible by physical arguments but could not establish them rigorously. 
It is the aim of the present paper to prove rigorously the existence of smooth solutions 
of problem (jl.ip . 

A problem closely related to (jl.ip has been treated in |16j . It is inspired by the 
dynamo problem with plane symmetry (which means d = 2) and differs from (jl.ip in 
that c = and by a different asymptotic condition at spatial infinity. The authors 
treat this problem and two related ones in arbitrary dimensions and prove existence of 
solutions and, moreover, some decay results. However, this problem does not precisely 
meet the requirements of the non-radial problem (jl.Sp nor those of the axisymmet- 
ric problem. Moreover, only weak solutions are established, whose behaviour at the 
boundary remains open. 

The basic idea of our treatment is to consider (jl.ip as parabolic problem in a 
bounded domain with non-local boundary condition, and to carry over the well-esta- 
blished methods for linear parabolic equations with standard boundary conditions such 
as Dirichlet's or Neumann's boundary condition to our situation (cf. e.g. [7]). In section 
2 we solve weakly a Poisson-type problem related to (jl.ip in all space. The regularity 
of the weak solution follows with standard arguments in G and G, only at dG we need 
special considerations. In section 3 we treat the corresponding eigenvalue problem, 
introduce the associated Fourier series, and characterize the elements of various function 
spaces which are useful in the remainder of the paper by the behaviour of their Fourier 
coefficients. In section 4 problem (jl.ip is solved by a Galerkin procedure and, for 
sufficiently smooth data and suitable compatibility conditions, the smoothness of the 
obtained solution is established. In two appendices the relation between problem (jl.ip 
and the non-radial-flow as well as the axisymmetric problem is elucidated. A third 
appendix collects some facts about exterior harmonic functions, and clarifies the relation 
between different kinds of spatial decay conditions - a topic about which there has been 
some debate in the literature. Finally, in a fourth appendix the completeness of the 
so-called poloidal free decay modes is proved: a commonly believed property which to 
our knowledge, however, has never been proved. 

2 A Poisson problem 

We establish in this section smooth solutions of the following Poisson problem in all 
space with suitable right-hand side /: 



Au = f 



in G 




u and Vu continuous 




Au = 



in G 



for Ixl — )• oo. 



m 



d 



(2.1b) 
(2.1c) 
(2.1d) 
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To obtain a weak formulation of problem ()2.ip let us multiply (|2.1b ) and ()2.1b ) with 
a test function v G C^{W^) and integrate over G and Br \G, G C Br, respectively. 
Assuming dG G one obtains after integration by parts: 

/ {Vu-Vv- fv)dx- / n-Vuvds = ^, (2.2a) 
Jg JdG 

/ _Vn-Vvdx+ / n-Vuvds- j^-Vuvds = 0. (2.2b) 

Here denotes a sphere with radius R and n is the exterior unit normal at dG. 
Adding up (|2.2b .) and (|2.2b ). one finds in the limit i? — )• oo that 



Vu-Vvdx= / fvdx (2.3) 



G 



In view of (|2.3p it is reasonable to consider functions satisfying the "finite energy 
condition" (cf. remark A.l) 



|Vt;|^dx =: < oo. (2.4) 

Condition (j2.4p together with condition (j2.1d ) motivate the definition of the real Hilbert 
space 

Ho :=clos{Co^(M'^), II -ll^} 

with scalar product {v,w)'}{ := J^d Vf • Vwdx. Observe that v G Hq is locally square- 
integrable. In fact, the Gagliardo-Nirenberg-Sobolev-inequality (cf. [71 p263]) implies 

„ „ ^„ „ 2d , , 

\Mlv(m.'1)<C\\v\\'h, p=j—^ (2.5) 



for any v E and a constant C depending only on d. Combining (|2.5p with Holder's 
inequality yields for any K ^W^: 

hh^K) < |i^|'/'||t^||Lni^) < c\K\'/%\\n ■ (2.6) 

So, defining (for later use) 

n:={ve HUR'') : ||t;||^ < oo}, 
there holds clearly Hq C H. Inequality ()2.6p implies, in particular, 

\\v\\l2(^g) < CgMIh , v^Ho (2.7) 

with Cg := C\G\^/'^- thus, v £ Tio yields in eq. (H^]) with / £ L'^{G) a finite right-hand 
side as well. 

A function u G TIq is now called a weak solution of problem (j2.ip with / G L'^iG) 
iff ()2.3p holds for all v eHo. Rewriting ()2.3p in the form 

(n, v)n = if, v)l2i^g) for any v ^Hq (2.8) 
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and noting that (/, •)l^(^g) defines a bounded linear functional on T-Iq (due to ()2.7p ) the 
existence of a unique weak solution follows immediately from the Riesz representation 
theorem. 

Concerning regularity of the weak solution let us define / : M'^ — t- M by / = / on 
G and / = on G. It is a standard result about interior regularity (cf. e.g. O p309f]) 
that / G H''{M.'^) implies u G i?f+^(M'='), k G Nq. This result means in particular: 

u G Hl^iR''), (2.9a) 
fGH\G)^uGHf^+\G), feGNo, (2.9b) 
u G C°°(G). (2.9c) 

With the usual Sobolev embeddings <^M )) implies u G C^{G) if / G H''{G), k > d/2. 
So, choosing suitable test functions in ()2.3p we find u satisfying ()2.1b ) in G and, for 
sufficiently regular /, ()2.1k ) in G. At dG, however, we need some finer considerations. 

As usual we flatten dG locally by means of a diffeomorphism $:[/—)• PV^, x i—)- y 
with $(xo) = 0, xo G U<r^^G, detZ)$ = 1 (cf. e.g. [f, p626f]). Setting v{y) := u{<l>~''{y)) 
and g{y) := /(^>~^(y)), <I> transforms eqs. (|2.ip localized on U into 

-Lv = g mW~, (2.10a) 

Lt> = in (2.10b) 

f and Vv continuous in W. (2.10c) 



Here, 



'yj 



is a uniformly elliptic operator with coefficients determined by $, W := W r\ {yd < 
0} = $(C/nG), Ty+ := Wr\{yd > 0} = $(C/nG), and := Wr\{yd = 0} = ^'(C/naG). 
A function v G -ff^(VF) is called a weak solution of problem ()2.10p with g G L?'{W~) iff 



gwdy 
iw- 



for any w G i/Q(VF). Defining "non- isotropic" Sobolev spaces H^^^{W), /c G Nq by 
Hj^^^{W) ■.= {v :W D'^v G L'^{W) for any multiindex a with \a\ <k,ad = 0} 



with norm 



\a\<k 



=0 



the regularity of v over 14^'^ is characterized by 

Lemma 2.1 Lei g G -?//^„„(M/^~) wzi/i k > {d+ l)/2, Zei v he a weak solution of ()2.10p . 
and /ei r 6e so sma/Z i/iai i/ie cylinder V := {y = {y' ,yd) G M*^ : < r, < 2r} d W. 
Then v G G^(Vi/2) with V1/2 := V Ci {\yd\ < r} and we have the estimate 

ll^llci(yi/2) - ^^^^^^^Htjw-) + hh^iw)) (2.11) 
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with C depending on W , V , k, and d. 
Proof: We prove first the estimate 

\\<HVy.) ^ ^ E E W^LD-vWl^^y^ (2.12) 

j=0,l,2 |a|<fc 

for functions v G C^~^'^{V) by applying Sobolev embeddings separately for y' and y^. 
Let us start with a 1-dimensional estimate on the interval (— r, r). For functions / G 
C"'^((— 2r, 2r)) we have 

|/(y)| < |/(y) - f{x)\ dx + ^ \f{x)\ dx 

|/'(x)|dx+- / |/(x)|dx. 

-2r ^'^ J-2t 

This implies that 

inax|/|2<C7( C |/'(x)|2dx+ C \f{x)\^d.x\ (2.13) 

for some constant C = C(r). On the other hand, Sobolev embeddings imply for 
k > {d+l)/2 

\\f\\l.0^<C\\f\\l.(^^,) = CY^ I l^"/(y')Pdy' (2.14) 

for some constant C = C{k,d,r) and B'^. := {y' G M'^"^ : \y'\ < r}. Combining ()2.13p 
and ()2.14p we obtain 

max ||.(,y.)||^,^-,^ ^ ^ E E r / |9;,^"^(y',y.)l^dy'dy. 

(2.15) 

i=0,l |a|<fc 

and analogously 

max ||5,,^;(.,y.)||^„(^,^ = E E H^^.^^Hli^cv)- 

' j=l,2 |a|<fc 

This is ([21^ . Now let v £ L^{V) with bounded right-hand side in ([232]) . It follows 
with standard arguments that v has a representative in C^{Vi/2)- 

Next, let g be the trivial extension of g on W. g G H^^j-^{W~) implies then g G 
^tani^)^ i.e. G L'^{W) for |q| < A;, = 0. From interior regularity for weak 
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solutions of ([211]]) follows D^^v G L'^{V) for < /c + 2, /3rf < 2 and the estimate 
j \D^v?<iy<c(Y^ j \D-9?<iy+\\v\\l.^^^ 

with C depending on W , V, k, and d. Together with ()2.12p this proves (j2.1ip . □ 

Remark 2.2 Observe for later use that lemma 2.1 applies in particular to the case 
g := Av|iy- with A G M. In fact, given any A; G N we find v\y[/- G Hf^j^{W~) for 
some V (S: W ^ W just by iterating the interior regularity argument. So, lemma 2.1 
applies to this eigenvalue problem as well and yields C^-smoothness of solutions over 
the boundary. 

Concerning the original problem ()2.ip . lemma 2.1 implies for dG G C^^"^ (i.e. ^> G C'^"''^) 
and / G H''{G), k> {d + l)/2 that u e over dG. Collecting the foregoing results 
we have 

Theorem 2.3 (Solution of the Poisson problem) Let G C M*^, d> 3 be a bounded 
domain with boundary dG and f G L^{G). The Poisson problem (j2.ip has then a unique 
weak (in the sense of eq. (|2.3p ) solution u G H^^^iW^). Moreover, if dG G G^^"^ and 
f G H''{G) withk > {d+l)/2, thenu is a classical solution, i.e. u G G'^ (R'^)nG^ (GUG) , 
satisfying pointwise eqs. ()2.ip . 



Remark 2.4 The condition d > 3 is crucial for theorem 2.3. In lower dimensions prob- 
lem (|2.ip is over deter mined; to obtain nontrivial solutions one has to relax condition 
()2.1b ) or (|2.1d ). As an example in M consider the problem —u" = 1 in G := (—1, 1), 
u" = in G, and lim|3,|_!.oo u{x) = 0. The only continuous solution is u{x) = ^ — ^x'^ 
in G and u = in G, which is not G^ over dG. Similarly, in M^, the only continuous 
solution of —An = 1 in G := {j; G : |x| < 1}, Au = in G, lim.\^\^^u{x) = is 
u{x) = I — ;j|xp in G and w = in G. 

On the other side, G^-solutions of these problems exhibit linear (in M) or logarithmic 
(in M?) growth at infinity, d = 3 is the lowest dimension, in which G^'^-solutions vanish 
at infinity. In fact, the analogous problem in has the G^-solution u(x) 
in G and u{x) = in G. 



1 1U|2 



2 



X 



Remark 2.5 Concerning the asymptotic behaviour of u for large x we refer to well- 
known facts about exterior harmonic functions (cf. appendix C). In particular, condition 
()2.1d ) implies the representation (jC.ip . which yields the asymptotics 

u = Yo\x\'^-'^ + 0{\x\^^'^), Yo^R for|x|^oo. (2.17) 

The decay is faster (viz. 0{\x\^~'^)), if / has vanishing mean over G. In fact, one obtains 
from ()2.ip . ()2.17p . and Gauss's theorem: 

\G\J = - [ Audx = - [ n-Vuds = - [ -(^ -Vuds = {d - 2)\Si\Yo. (2.18) 

Jg JdG Jsr m 
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Here, / := J^, / dx and 5"^ C G; the last equality arises in the limit i? — oo. So, 
/ = implies Yq = 0. 

Let us note as an aside that in general n = can only be achieved at the expense 
of relaxing (|2.1d ) to u{x) = 0(1) for |x| — )• oo. In this case the function u — u with u 
being a solution of (|2.ip is obviously a zero-mean solution ( see also remark 2.6). 

Remark 2.6 If G is a ball Bji it makes sense to consider spherically symmetric so- 
lutions of problem (|2.ip . i.e. {u){r) = for any r > 0. Obviously, {u) = implies 
(/) = 0; on the other hand, any solution of ()2.ip with spherically symmetric / is spher- 
ically symmetric. This follows from the unique solvability of the sub-problem for the 
spherical mean arising from (IZTdI^ 

-{{u)"+id-l)/r{uy)={f), 0<r<R, 
{u){R) = Yo/R''-^ (n)'(i?) = (2 - d)Yo/R''-\ 

Observe here the representation (|C.ip . which implies (u) = YQr'^~'^ for r > R, and 
(|2J8]1 . which implies Yq = if (/) = 0. 

Note, finally, that spherically symmetric solutions decay at least like for 
large x. 



3 The Eigenvalue problem 

We treat in this section the eigenvalue problem corresponding to the Poisson problem 
(|2.ip of the previous section: 

—Au = Xu in G, (3.1a) 

Au = in G, (3.1b) 

u and Vii continuous in M'^, (3.1c) 

u{x) — )• for |x| — )• oo. (3. Id) 

According to theorem 2.3 we have for any / G L'^{G) a unique weak solution u G 
of (|2.ip , defining thus a Green operator 

g:L\G)^no, f^u. 

Using dST]), dSSD we obtain 

116^/11?^ = {u,u)n = (/,^)l2 < WfWL^Qfh^ < CGWgfWnWfh^ (3.2) 

and 

(/,^/)l2 = {f,u)L2 = {u,u)n = WGffn > 0. 

'^Prime means differentiation with respect to r. 

■^The symbols and without specified domain mean always L^{G) and H''{G), respectively. 
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Therefore, ^ is a bounded linear operator between Hilbert spaces, which is, furthermore, 
positive and hence symmetric. Restricting u on G one obtains the operator 

g:L^{G)^L\G), /^uIg, 

which is hkewise bounded and symmetric, and, moreover, compact due to the Rehich- 
Kondrachov theorem and the observation {u\g '■ u G T-Iq} = {u\g ■ u £ Ti} = H^{G). 
The spectral theorem for symmetric compact operators in Hilbert spaces establishes 
now a complete (in L?'{G)) orthonormal system {u„ : n S N} of eigenvectors of Q: 

gvn = X-\n, nGN, (3.3) 

with real, positive eigenvalues A^^ of finite multiplicity and lim„_s.oo = 0. In order 
to solve the original problem let us define the "harmonic extension" Vn ■= ^nQvn of 
Vn on M'^jfl By definition the Vn are weak solutions of the Poisson problem (|2.ip with 
/ := Vn and u = X~^Vn': thus, eq. (|2.8p takes now the form 

(vn,v)n = A„(w„, v)l2 for any v £ Tio , (3.4) 

i.e. the pair {vn,Xn) is a weak solution of the eigenvalue problem ()3.ip . Note, finally, 
the uniqueness of the harmonic extensions in T-Lq, which is implied by the uniqueness 
of weak solutions. 

Concerning regularity we have Vn G G°°{GuG)r\G^(M.'^), which follows from ()2.9b ). 
iterating ()2.9b ). and remark 2.2. So, (vn, An) is also a classical solution of problem ()3.ip . 
We summarize these results in 

Theorem 3.1 (Solution of the eigenvalue problem) Let G C W^, d > 3 be a 

bounded domain with -boundary, k > [d + 5)/2. The eigenvalue problem ()3.ip then 
has a countable set of eigensolutions {{vn,Xn) : n G N} satisfying ()3.4p . and their 
restrictions {vn : n G N} constitute an orthonormal basis of L?'{G). 

Powers of the inverse Green operator and their domains of definition turn out to 
provide the right setting for the solution of the evolution problem in Section 4. The 
elements of these spaces can be characterized by the decay behaviour of their Fourier 
coefficients when expanded in the above eigenfunctions. This motivates the 

Definition 3.2 Let {vn : n G N} be the complete orthonormal system with associated 
eigenvalues A„ according to theorem 3.1, and a G M. We define then the space of 
"formal series" 



S := < ^CnVn : c„ G M L 

n=l 



with non-negative functional 



oo ^ oo \ 

:cS^[0,Oo], Y^CnVn^ lY.\l''\Cnn 
n=l ^n=l ^ 



1/2 



^In this section a quantity with tilde always means the "harmonic extension" on R'* of a quantity 
defined on G. 
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linear mapping 

oo oo 
n=l n=l 

and subspaces 

Da := -D(^") = e 5 : < oo} C 5. 

Obviously, A" maps L»„ into -L2(G). Further more, there is Dq, C D/j if a > /? and 
L>o = L'^{G). Thus, if a > 0, there is Da C L'^{G) and v G has the representation 

oo 

V = ^^{Vn-,v)L2Vn. (3.5) 

n=l 

With the pairing 

oo oo oo 

{w, v) := ^ dnCn for tt; = ^ d„v„ G w = ^ c„7;„ e (3.6) 

n=l n=l n=l 

D.Q, is the dual space of Dq,. 

Applying A on (finite) linear combinations of f„ we find ^ = G"^, and -D(^) turns 
out to be the maximal domain of definition of Similarly, -D1/2 is related to T-Lq. 

More precisely, we have 

Theorem 3.3 Let G C , d > 3 be a bounded domain with C°° -boundary dG and 
let {vn : n G N} be the complete orthonormal system defined by the eigenvalue problem 
()3.ip . Let, furthermore, A and Da be as defined in definition 3.2, and Q be the Green 
operator associated to the Poisson problem ()2.ip . Then, 

Do = H^{G) = L\G), 

Di/2 = {v\g ■ V G Ho and v\q is harmonic} = H^{G), (3-7) 
i.e., in particular, any v £ -D1/2 has a unique harmonic extension v G Tio, and 

D^ = G{L\G)) = {ve H\G) : v G HI^{W'')}. (3.8) 

Higher order spaces are characterized by 

Dk/2 = e H'^iG) : A^v G HI^{R'') for i = I, . . . , [k/2]} A: G N \ {1} , (3.9) 

where w G Tio again denotes the harmonic extension of a function if G 1)1/2 = H^{G) 
and [r] := max{j G N : j < r} is the integer part of r. 
On Df./2 we have the equivalence of norms: 

II • IU/2 ~ II • Ib^ , /cGNo. (3.10) 

To prove the theorem the following lemma is helpful. It improves the regularity result 
()2.9b ) and provides the pertinent estimate. 
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Lemma 3.4 Let G C M*^, d > 3 be a bounded domain with C''^'^ -boundary dG and 
f G H^{G), G Nq. Let, furthermore, u G 1-Lq be the weak solution of problem ()2.ip . 
Then u G H^^'^{G) and we have the bound 

< = CIIA^zll^.(G) (3.11) 

with a constant C depending on G, k, and d. 

Proof: The case A; = is already implied by the interior regularity result ()2.9b ). In 
fact, u G Hf^^{W^) means (see 0, p309) 



|h2(g) < c{\\f\\mK) + \\u\\mK)) , (3.12) 



where / denotes the trivial extension of / onto M'^ and K some bounded domain such 
that G K. Combining (|2.6p with p.2p we obtain 



\\u\\l2i^k) < CK\\u\\n < Cii-Cc 11/ II L2(G) > (3-13) 
and thus ()3.12p takes the form 

MH^(G)<C\\fh2iG)- (3.14) 

No boundary regularity is required for this result. 

The case k > needs separate considerations of tangential and normal derivatives 
at dG. We refer in the following to the situation, where dG has already been flattened 
as explained in the paragraph before lemma 2.1 and we use the notation introduced 
there. So, given g G L'^(W~) we assume v G H^{W) to be a (weak) solution of 

/ cLijdy^v dy^w dy = / gwdy (3.15) 
ij=iJw Jw- 

for any w G Hq{W). Let g be again the trivial extension of g onto W. Now we assume 
higher tangential regularity of g, i.e. D'^g G L'^{W~) for |a| < k, ad = 0, which implies 
D°^g G L'^{W). From interior regularity for weak solutions it follows that D^^v G L'^(V) 
for \I3\ < k + 2, fSd < 2 and any V W, together with the estimate 

[ \D^vfdy<c(Y, [ \D^g\'dy+\\v\\l2^^X (3.16) 
As to normal derivatives, note that higher interior regularity implies 

d 

-D" j;a,,(a,,a,^,t;) = D"5 (3.17) 
to hold a.e. in W~ . Writing ()3.17p with a = (0, . . . , 0, 1) in the form 

d 

adddy^v = - dy^dy^{aij dy^v) - 2dy^adddy^v - d^^adddy^v - dy^g , (3.18) 

i+j<2d 
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we find, by uniform ellipticity, dy^v to be bounded in by the right-hand side of 
(|3.18p . which is at most of second order in dy^v. So, ()3.16p may be apphed and we 
arrive at 

l^_\dlv\'dy<c( Yl l^_\D"9\'dy + \\v\\l,^^X (3.19) 

|a|<fc 

The case of arbitrary higher derivatives is now easily proved by induction. So, we find, 
finally, that p.l6p holds without restriction on ad and /3d, respectively. 

To complete the proof one has, as usual, to cancel the change of variables, to cover 
G by local patches, to sum up the corresponding local estimates, and to use once more 
(I3J3]1 . □ 

Proof [of theorem 3.3]: The case A; = is trivial. To prove ()3.7p let v G D1/2 be 
decomposed as in p.Sp . i.e., v = Yy'^=ii'^n,v)L2 Vn with ^^"^=1 ^n\{vn,v)L2\'^ < 00. We 
define 



V := 'Y{Vn,v)L2 Vn 
n=l 

with Vn G T-Lq being the unique harmonic extension of Vn- Computing 

00 00 

n,m=l n=l 

where we used (|3.4p . we find v G Hq. In order to prove that v\q is harmonic it suffices 
to show {v,v)-}i = for any v G C^{G) with suppu C G; this follows immediately with 

00 00 

{V,v)n = Y{Vn,v)L2{Vn,v)n = Xn{Vn, v) {Vn,v) ^2 = 0, 
n=l n=l 

since {vn,v)i2 = J^VnVdx = 0. Therefore, v G {v\c '■ v G Hq and v\q is harmonic}. 
The opposite inclusion follows again with (j3.20p . 

The inclusion {v\c ■ v £ Hq and vjg is harmonic} C H^{G) is obvious; the opposite 
inclusion follows with theorem C.l: Let v G H^{G), vq an //^-extension of v on with 
bounded support, w := vq\q G H and u the exterior harmonic solution of ()C.6p . The 
function u defined by u in G and n in G is then the sought-after function G Tio- This 
proves ()3.7p . 

To estimate the 1/2- norm we supply the above construction with bounds: Let 
supp Vq C K, then 

\\vo\\m{Rd) < G{G,K)\\v\\hi(^g) ■ 
Thus, using the minimizing property of solutions of ()C.6p . 

\\u\\^ < \\w\\^ < \\vo\\n < G{G,K)\\v\\hi(g) > 

and with (iOOjl : 

II ||2 ll~l|2 _ IIY7 ||2 I ||^||2 ^ r^w ||2 

lPlli/2 - W^Wh - II^'"IIl2{g) + - '-'W'^WmiG) ' 
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which is one half of ()3.10p for k = 1. The other half follows by ()2.7p . 

We show next that Di = GiL"^). Let u € g{L'^) and / = E^=iK, /)l2 ^-n e L"^ 
such that u = Qf. Computing the coefficients of u we find with (jS.Sp : 

oo 
n=l 

Thus, 

oo oo 

m=l m=l 

hence u & Di. If, on the other hand. 



'^{vn,u)L2Vn and A^|(fn,n)^2p < oo, 

then 



u = 

n=l n=l 



n=l 



is well-defined and we find Qf = u. 

The inclusion Q{L'^{G)) C {f G H^{G) : v G Hf^^{M.'^)} is an immediate consequence 
of the if ^-regularity of weak solutions. To prove the opposite inclusion let w G H^{G) 
with harmonic extension w G %q H Hf^^{W^). Defining / := —Aw G the Poisson 
problem (j2.ip yields a solution S G ?^on-ff;Qg(M'^). So, we have pointwise a.e. /S.{w—u) = 
Q m.W^ ioT w -u £ ^-LQ<r^ H^^iR'^), which means w — u is harmonic in (by Weyl's 
lemma) and, moreover, w — u = (by Liouville's theorem). Thus, w = u and, in 
particular, w = u = G{f)- 

To estimate the 1-norm of u G D{A) observe that for its harmonic extension holds: 
V ^%Qr\ Hf^^{W^), and for the eigenfunctions Vn- Vn G C^(M'^). So, by (f3^ we can 
calculate 

- {XnVn,v)L2(G) = - Vv„ • Vvdx = / Vn Av dx = {Vn, Av) l2^^g) (3-21) 

and therefore obtain 

oo oo 
\i = =Y.^n\iVn,v)i^2\' = X \{Vn,Av)L2\' = \\Av\\l, , (3.22) 



n=l n=l 



which implies ||f ||i < C||'y||/i'2(Q) with a constant C depending only on d. To prove the 
opposite inequality we combine ()3.22p with ()3.1ip t— n: 

llt'lli = IIA^II > ^\\v\\h^g)- 

This proves ()3.10P t.-9. 

The case A; > 2 is proved by induction. Let v G -0^/2+1, A; G N, i.e. Av G By 
assumption we have Av G H^{G) and 

/^^v G Hl^iR'^) for i = 1, . . . , [A;/2] . (3.23) 
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(Note that for v £ condition ()3.23p does not yet make sense and can be omitted.) 
By ()3.22p the condition Av G H^{G) means A?; G H^{G), and lemma 3.4 imphes 
V G H^^'^{G). Moreover, we have v G Hf^^iW^), which complements condition (|3.23p . 
So, we conclude 

^; G {z; G H^-^^{G) : A^^v G Hl^iW^) for i = 1, . . . , [k/2] + l}. (3.24) 

To prove the opposite inclusion let v be as in p.24p . We set to := Av and have by 
assumption 

w G |v G H''{G) : G HIcO^"^) for i = 1, . . . , [fc/2]} = Dk/2 ■ 

Computing the k/2 + 1-norm of v we find with (|3.2ip 

oo oo 

ll^'llfc/2+1 = ^>^n\>'n{vn,v)L2\'^ = ^ (-?;„, u;)i2 ^ = ||u;||^/2 < (3-25) 

n=l n=l 

and thus, v G -Dfc/2+i- This completes the proof of ()3.9p . 

As to the equivalence ()3.10p we proceed likewise by induction. Assuming v G 
Dk/2+i-, A; G N we find by ()3.25p and by assumption 

ll^^llfc/2+l = l|Al'||fc/2 < C\\/\v\\Hk < C\\v\\Hk+2, 

whereas the opposite inequality follows by (|3.1ip : 

ll^^llfffc+a < G\\Av\\h^ < C\\Av\\k/2 = C'||l'||fc/2+i. 

This completes the proof. □ 

Remark 3.5 Iterating ()3.2ip one finds on for integer values a the following alter- 
native formulation of the a-norm: 

\\v\\k = \\A''v\\L2^G), wG-Dfc, k£N, 
and by (|3.2Up for half-integer values: 

ll^^llfc+i/2 = l|VA'=u||^2(Rd) , veDk+i/2, fcGNo. 

4 The evolution problem 

We solve in this section the evolution problem (jl.ip by means of the spaces provided 
in the last section. According to eq. ()3.7p . v G -^1/2 has a harmonic extension v on 
M*^; so, when working with these spaces it is sufficient to consider problem (jl.ip on the 
simpler domain G x However, a nontrivial asymptotic function u^o does not fit 



14 



into this framework. Therefore, in a first step, is ehminated by the time-dependent 
shift u — Uoo '■= Us- In terms of Us problem (jl.ip reads 



dtUg — aAug = b ■ Vug + c Ug + f 


in G X M+, 




(4.1a) 


Aus = 


in G X R+, 




(4.1b) 


Us and Vus continuous 


in M'^ X M+, 




(4.1c) 


Us{x, t) — )• 


for X — )• oo, 


, i G M+, 


(4.1d) 




on G X {t = 


0} 


(4.1e) 



with f := cuoo - ^Uoo- 

The "simphfied" problem takes then the form 

v = -aAv + Bv + f, (4.2a) 
^^(0) = vo, (4.2b) 

with the operator A as defined in Definition 3.2 and the lower-order operator B defined 
hy Bv := b- Vv + cv. Here v is a mapping from [0, T), T > into some function space 
over G. As explained above a reasonable such space is = H^{G) with at least 

k = 1 (cf. theorem 3.3). Moreover, for T > 0, when starting with vq G ^1/2 and taking 
into account parabolic smoothing we expect v G L^((0, T), Di) which means in view of 
(lOk ) i G L2((o,T),L2(G)). This motivates the 

Definition 4.1 Let T > and vq G ^ function v G L^((0,T),Di) iwii/i weak 

time derivative v G L'^{{0,T), L'^{G))) satisfying (jOk ) as equality in L'^{{0,T), L'^{G)) 
and (|4.2b ) as equality in -D1/2 called weak solution of problem (j4.2p . 

Condition ()4.2b ) makes sense for weak solutions due to the following interpolation 
result: 

Lemma 4.2 Lei G be a bounded domain with smooth boundary, T > 0, and k G 'Nq. 
Let, furthermore, v G L'^{{0,T), H''+^{G)) and v G L^{{0,T), H''-^{G)). Then 

vGC{[0,T],H''iG)); 

moreover, the mapping 1 1— )• ||^(^)||/,2(c') is absolutely continuous with derivative 

l^^HmUG) = {m^<t)) (4-3) 

for a.e.te [0,T]. 

For a proof we refer to [3 p287f]. We note only in the case k = that -D_i/2 is the 
dual space of -D1/2 = L[^{G), thus -D_i/2 C H^^{G). {■ , •) denotes the dual pairing as 
defined by <\3.6\i . 

Theorem 4.3 (Weak solution of the evolution problem) Let T > and a,b,c E 
G(Gx [0,T]), a>ao> 0. Let, furthermore, vo G -D1/2 and f G C{[0,T], L'^{G)). Then 
problem (j4.2p has a unique weak solution v. 
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Proof: We start with the construction of Galerkin approximations using the complete 
system {v,^ : v G N} from theorem 3.1. Let Py be the orthogonal projection in L'^{G) 
onto span{^;,,}, := 0"=i P^., and let G be the unique solution 

of the following finite-dimensional initial-value problem 

^^(n) ^ p{n) ^ ^(n) ^ ^ ^(n) ^ y.) ^ (4 4^^ 

i;(")(0) = pWt;o. (4.4b) 

Note that P^") commutes with ^ but not with a or B. From standard results about 
ordinary differential equations follows v^^^ G C^([0, T], Di) for any n G N. 

Next we derive some a-priori estimates for f^"^ uniform in n, which allow to extract 
a weakly convergent subsequence of the sequence (v*-"^) of Galerkin approximations. We 
first show that 

max||t;(")||2/2<C = C[t;o,/;r]. (4.5) 
Taking the scalar product of (|4.4b ) with Av^'""^ we obtain 



< -ao|b(")||? + (^?;("\Sr;("))L2 + (^?;("),/)i2. 

Observing that ;B is a bounded operator from C([0, T], D1/2) into C([0, T], L2(G)) there 
is a constant Ci such that 

(^^;W,5t;W)i2 < (Ciao)l/2||^(n)||^||^(n)||^^^ < ^ ||^(n)||2 ^ ^ • 

Setting max[o^T] Il/lli2 =: C2ao we thus obtain 

^^\\v^''^\l/2<Ci\\v^^X/2 + C2, (4.7) 
and Gronwall's inequality yields 

\\v^-X/2 < e^^*ll^(")(0)||?/, + |(e^^* - 1) < e^^^ll.r II?/, + |(e^^^ " l) 

on [0,T], and hence ()4.5p . 

To obtain a bound on v^""^ in L^((0, T), Di) we estimate similarly to (|4.6p : 



\ i < -«olb(")||? + ^lb(")||? + Ci||.W||?/, + |||.W||? + . (4.8) 

Using ()4.5p we rewrite ()4.8p in the form 



aolb(")||? < \\v'^''^\\l/2 + 2(CiC + C2) ; 
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thus, integrating over [0, T] and observing ()4.5p once more yields the bound 

r \W''^\\ldt<C . (4.9) 

^0 

With (gJD the right-hand side in (Ob ) is obviously bounded in L'^{{<d,T),L'^{G)), i.e. 
there is C such that ^ 

||^^i(")||'2dt<C (4.10) 

for the sequence of (classical) derivatives (;^^^^"^)- 

The bounds ()4.9p and ()4.10p imply that there is a subsequence {ni) and functions 

V e L2((0,r),Di), {; G L2((o,r),L2(G)) such that v is the weak limit of (u^"')) in 
L2((0,r),Di) and v that of (^w^"'^) in ^^((0, T), L2(G)), respectively, and moreover 

V is the (weak) derivative of v. 

Testing (lOh ) with functions w of the form ty(t) = X;™=i f^i.(i) i";. e ^^([0, T],L'^{G)), 
where m <n and djy : [0, T] — )• M are smooth functions, and integrating over [0, T] yields 



Jo "i* Jo 

rT 

= / (-a^?;(") +/,-«;) 2 dt. 
Jo 

Setting n = rii we find in the limit / — )• oo that 

/ {v,w)L2dt= [ {-aAv + Bv + f,w)L2dt. (4.12) 
Jo Jo 

Since test functions of this type are dense in L^((0, T), L^(G)), eq. (|4.12p holds for any 
w e L2((o,T),L2(G)). This proves g2Ji) to be an equality in L'^{{0,T), L'^{G)). 

Inserting w G C^{[0,T], L'^{G)) with w{T) = in ()4.12p . after integration by parts 
on the left-hand side we find 

Doing the same in (|4.1ip yields in the limit n/ — )■ oo on the left-hand side, 



(4.11) 



T 



d 



L2 



Since w{0) G Lp'{G) is arbitrary we have u(0) = vq in L^{G). These results prove v to 
be a weak solution of problem (|4.2p . 

Finally, to prove uniqueness of the weak solution consider vi — V2 =■ vq satisfying 

vo = -aAvo + Bvo , vo{0) = . 

Setting in eq. v := vq e L^{{0,T), D{A)), w := Avo, and / = we obtain for 

a.e. t G [0,r]: 

{vo,Avo)l'2 = -iaAvo,Avo)L2 + {Bvo,Avo)l2 ■ (4.13) 
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The left-hand side of ()4.13p combined with ()4.3p takes the form 

{vo,Avo)l'^ = {A^/'^vo,A^/^vo) = ^ ^ WA^^^oWh = ^ ^ \\vo\\l/2 > 

whereas estimates analogous to that leading to (|4.7p show that the right-hand side of 
(|4.13p can be bounded by ^Ci||wo||^/2- Since 1 1— )• ||fo(i)||i/2 is absolutely continuous by 
lemma 4.2, applying Gronwall to the inequality 

^ \\vo\\l/2 < Cilko||?/2 for a.e. t £ [0,T] 

with 1 1 I'D (0)1 1 1/2 = yields the desired result vq = 0. □ 

Remark 4.4 As to the original (shifted) problem ()4.ip theorems 3.3 and 4.3 imply 
that Us{x,t) := \v{t)]{x), where v{t) is the harmonic extension of v{t), satisfies ()4.1b .b) 
for a.e. {x,t) € R'^ x [0,T]. 

Higher regularity of the weak solution depends on the smoothness of the coefficients 
and the initial-value, and suitable compatibility conditions among these data. There 
holds 

Theorem 4.5 (Higher regularity) Let T > 0, e N \ {!}, and a,b,c e Cf^(G x 
Let, furthermore, vq £ -D(a..+i)/2; -a{-,0)Avo + B\t=oVo + f{0) G -D(fc_i)/2, and 
f G C^{[0,T], H^{G)). Then the weak solution v of problem (j4.2p satisfies 

veL\{0,T),D,/2 + ,), veL\{0,T),D,/^), v e L\{0,T), D,/^^,) . 

Proof: Higher spatial regularity is easily obtained via the operator A: applying A^^"^ 
on ()4.2p and setting A''^'^v=:w, 

j^k/2j ^. ^ik)^ A^'^BA-^''' =: B^^\ (a - A^'''aA-^''')A =: 

and A^^'^vq := wq, we obtain 

w = -aAw + B^^'^w + 6^'^w + f^^\ (4.14a) 
w{<d)=wo. (4.14b) 

/'^'^^ and wq fulfill the prerequisites of theorem 4.3, 

i3(^-):C([0,r],Di/2)^C([0,r],L2(G)) 

is again a bounded operator, and C^^^ is of the same type as B'^^\ Thus theorem 4.3 
applies to (|4.14p with the result 

w G l2((0, T), Di) , we L\{0, T),L\G)) , 



''in this notation the upper index at "C" refers to the order of spatial derivatives and the lower one 
(omitted if zero) to temporal derivatives; so, a G Ci (G x [0, T\) means a, dta, and -D"a, |a| < k are all 
continuous functions on G x [0,r]. 
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I.e. 

V G L2((0,T),I)fc/2 + i) , V G L2((0,T),Z)fc/2) . 

To obtain higher temporal regularity we need some more a-priori estimates for the 
Galerkin approximations ti^"^ . Note that a, 6, c are at least G Cl{G x [0,T]), t;o G Dg/s, 
and / G C^([0, T],D{A)). So, inserting := Av^'^'^ into (|iHk l and taking the scalar 
product with Aw^'^^ we obtain 

This is analogous to ()4.6p and the subsequent estimates leading to ()4.5p yield now 

max 11? = max ||?;(") ll^ ,„ < , (4.15) 
[o,T] ^'"^ [o,r] •^'"^ 



and, via (|4.4b ). 



max||-^?;(")||^.^ < C4. (4.16) 
[o,r] "dt "V2 

On the other hand, differentiating (|4.4b .) with respect to t, setting 
we obtain 

^^(n) ^ p(n)( _ a^^W + + / _ dtaAv^''^ +Bv^''^) , (4.17) 

which is of type (j4.4b ). Complementing ()4.17p by the initial- value i>(")(0) = P^^^vq, 
where vq := —a{-,0)AvQ + B\t=oVo + /(O) G -Di/25 we obtain again a bound of type 
(|4.5p . now on Next, we modify for (|4.17p the argument which leads from (|4.5p to 
(|4.9p . Taking the scalar product of (|4.17p with Av^"^^ and observing that B is of the 
same type as B, thus using the bounds 

\\Bv^-^h2 < (Ciao)i/'||t'(")||i/2, \\Bv^^^h2 < (Ciao)i/2||^(")||i/2 

as well as max^QX] ^ C^, ()4.5p . and a > ao, \dta\ < A, max[o^T] 11/11^2 < ^2001 

we obtain 

I ^ 11^^"^ < -«o||^(") II? + |||^(") II? + 2 II^W ||?/2 + 1 1|^(") II? + 2 C2 

+ ^||^W||? + 2^C75 + ^||^(")||? + 2CiC. 
8 ao o 

This estimate is analogous to ()4.8p . so we have 

rT 



and, using (j4.17p once more. 



/ ||t'(")||?dt < Q 
Jo 
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Recalling the reasoning after ()4.10p there is thus a subsequence (n;^) =: (m) of (n) 
and a function v G L^((0, T), L^(G)) such that v is the weak limit of {■^jiv^'^^) and the 
weak time derivative of v. Finally, inspecting again ()4.17p we find that there is enough 
regularity of the data a, b, c, /, and vq left to improve the spatial regularity of v by 
the order k — 2, i.e. we have v £ L'^i{0, T),Df^/2 - i)- 1^ 

In view of lemma 4.2 and Sobolev's embedding theorems, theorem 4.5 implies the 
existence of smooth solutions. The following corollary formulates for the original evolu- 
tion problem (jl.ip sufficient (and not necessarily sharp) conditions in terms of classical 
derivatives for existence of classical solutions. In particular, to express the compatibility 
conditions in classical terms, note that u G C^{G) and u € C^(M'^) imply u G Hf^^{W^). 



Corollary 4.6 (Classical solution of the evolution problem) Let G C M°', d > 
3 be a bounded domain with C'^^^^^ -boundary, k > l + d/2, uq G C''^^{G), and a,b,c£ 
Cf(Gx [0,T]), Uoo G C^{[0,T]) for any T > 0. Let, furthermore, uo — Uoo{0), A^uq, and 
A*^^(aoAno + bo • Vuq + cqUq — Uoo{0)), i = I, ■ ■ ■ ,[{k — l)/2], where aQ = a(-, 0) etc., 
all -match to their harmonic extensions. Then problem (jl.ip has a unique classical 
solution u, i.e. u G Cl{G x M+) n C^{G x M+) satisfies pointwise eqs. (jl.ip . 

Proof: Fixing some T > 0, by theorem 3.3, lemma 4.2, and Sobolev's embedding 
theorems, from theorem 4.5 we find 

V G C{[0,T],H''+\G)) C C{[0,T],C\G)) , 

V G C{[0,T],H''-\G)) C C{[0,T],C(G)) . 

Thus, setting Us{x, t) := [v{t)]{x) with v{t) being the harmonic extension of v{t) G -Di/2) 
we have Ug G Cf{G x [0,T]) satisfying (|4.1k ) and ()4.1b ). Since v{t) G Ho is harmonic 
in G conditions (|4.1b ) and ()4.1d ) hold for Ug as well. To prove ()4.1b ) note that this 
condition holds for Vn- Recalling the maximum principle for harmonic functions this 
implies 

max I (x) I = max | (x)| 

for any multiindex a with \a\ < 1. Setting Smn{x, t) := J2u=m ^^^i^) Vu{x) we have then 
with Sobolev and relation (j3.10p 

max > max iL'^S'mnfx, t)| = max > max |L>"5m,„(x, t)| 

<C'^maX \\Smn{--,t)\\Hk(^G) ^ C\\Smn\\c{[0,T],Dk/2) ■ 

(4.18) 

So, convergence of v{t) = ^^1^=1 Cn(0 Vn in C([0, T], D;-/2) implies convergence oiv{t) = 
En=iCn{t)vn in C {[0,T], C^W")) , i.e. Us G C\R^ X [0,r]). 

Similarly, fixing any K (£ G and using the interior derivative estimate (cf. [9l p23]) 

iiiax\D^Vn{x)\ < C 'nia_x\vn{x)\ < Cmax |?^„(x)| 
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with |a| = 2 we find G C^iK x [0,r]) and hence Us G C^iG x [0,r]). 

Observing, finally, that T > is arbitrary we have in conclusion that Ug is a classical 
solution of problem ()4.ip and hence u := Ug + Uoo is a classical solution of problem (jl.ip . 
□ 

Remark 4.7 In d = 3 we may choose k = 3 and the compatibility conditions amount 
to 

uo ^^^(0) , A^To G C\m:^) (4.19) 

and 

(aoA'Uo + 6o-Vuo + cono-Uoo(0))~ gC^(m3). ^4 ^0) 

So, in the = admissible initial values uq are for instance C^(G)-functions 

with dl^uolgfj = 0, z = 0, . . . , 3, where 5„ denotes the normal derivative at dG. In the 
case Mo = "^^oo = const > 0, which was interesting in applications |llj . condition ()4.20p 
requires the coefficient cq to have a C^-smooth harmonic extension. 

Appendices 

A The non-radial-flow problem 

In the framework of magnetohydrodynamics the kinematic dynamo problem reads [1]: 



{r]V X B) = V X {v X B) , V -B = 


in G X M+, 




(A.la) 


V X B = 0, V .B = 


in G X M+, 




(A.lb) 


B continuous 


in X M+, 




(A.lc) 


B{x,t) = 0{\x\-^) 


for X —7- 00. 


, t G M4 


(A.ld) 


B{;0)=Bo 


on G X {t = 


0}. 


(A.le) 



Here, the induction equation (jA.lb ) describes the generation of the magnetic field B 
by the motion (with prescribed flow field v) of a conducting fluid (with conductivity 
> 0) in a bounded region G C M^. Outside the fluid region there are no further 
sources of magnetic field. Thus, B continues in G = M'^ \ G as a vacuum field and 
vanishes at spatial infinity. 

If G is a ball Br (or a spherical shell) the so-called poloidal-toroidal decomposition 
of solenoidal fields is especially useful [H [H] : 

B = Bp + Bt = -V xAS-AT, {S) = {T)=0. 

A denotes here the non-radial derivative operator A:=xxV,A-A=:£is the Laplace- 
Beltrami-operator on the unit sphere Si, and (•) denotes the spherical mean. The 
poloidal and toroidal scalars S and T, resp., are uniquely determined (e.g. in L?'{Bpi)) 
hyB: 

x-B = -CS, x-V X B = -CT. 
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In the following we refer to P := x ■ B instead of S as the poloidal scalar. 

In the case of a non-radial flow field, i.e. u • x = 0, and spherically symmetric 
conductivity, problem (jA.ip implies the scalar sub- problem p.3p for P: ()1.3h ,) is just the 
radial component of the first part of (jA.lb ) and ()1.3b ) is obtained by applying A on the 
first part of (jA.lb ). Condition (jl.Sb ) is in fact enough to ensure a continuous magnetic 
field B, i.e. continuous second-order derivatives of S" = —C~^P, since B involves not 
more than one radial derivative of S. The equivalence of (|A.ll i) with (|1.3d ) is clear for 
B ■ {x/\x\) and follows for the non-radial components with the divergence-constraint. 

Remark A.l In the mathematical treatment of problem (jl.ip it turns out to be useful 
to consider functions v satisfying the integral condition Jj^^j |Vw|^ dx < oo. In the 
context of problem (jl.Sp this condition can be interpreted as one guaranteeing finite 
total magnetic energy. In fact, the total energy of the poloidal magnetic field reads 

E[Bp] = \ f \V X AS\'^dx = l [ \VC^/^S\'^dx, 

2 Jm.3 2 Jr3 

and with the variational estimate 

inf "^^"^^(^^) = 2 
f^o,{f)=o ||/||l2(5^) 

one obtains the bound on E[Bp]: 

E[Bp]<-l \VCS\^dx = -l' iVPl^dx. 

B The axisymmetric problem 

The central assumption is here an axisymmetric magnetic field with representation 

B = VPx V4> + AV(p = --d^Pep + -Ae^ + -dpPe^ (B.l) 

P P P 

by two scalar quantities, the poloidal one P and the toroidal or azimuthal one A, 
depending (besides on t) on p and z with (p, (p, z) being cylindrical coordinates in 
\ {/^ — 0}- Inserting (|B.ip into the dynamo equation (|3.1ip the following sub- 
problem for the poloidal scalar P arises [3} 112^ [TO]: 

dtP - T]A^P = -vVP 
A*P = 
P and VP continuous 
Pip,z,t)^0 
\VP{p,z,t)\=0{p) 

-vp{p,z,t)\=o{{p'+z'r'/') 

Pi-, -,0) =Po, Po satisfying (B.2d,e) 



in G2 X R+ 




(B.2a) 


in G2 X M+ 


? 


(B.2b) 


in H X M+, 




(B.2c) 


for /5 — 7- 0, 


[z,t) eRxl 


(B.2d) 


for /5 — )• 0, 


>,t) G M X I 


(B.2e) 


for + 


00, t e M 


f, (B.2f) 


on G2 X {t 


= 0}. 


(B.2g) 
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A* is the elhptic operator A* := d'j, — ^dp + d1 on the half-plane i/ := M_|_ x M. G2 C H 
is the "cross-section" of some bounded region G3 C M^; more precisely, G3 \ dG^ with 
G3 := G2 X {0 < (p < 27r} is a bounded domain in M'^ with smooth boundary. Note that 
the axisymmetric flow field need not be solenoidal, the azimuthal component, however, 
w.l.o.g. can assumed to be zero. Condition (]B.2fe ) ensures a finite magnetic field on the 
symmetry axis {p = 0}. It implies the limit limp_^o ") ^) = Ps{t), where Pg depends 
only on t. As Pg does not affect the magnetic field it is set to zero for simplicity 
(condition ()B.2d )). Note that in [12l[10] conditions ()B.2d .e) are replaced by 

= 0(p2) for 0, G M X M+. (B.3) 

The cautious Backus [3J requires ()B.2b ) and (|B.3p . In fact, conditions ()B.2l i.e) imply 
(jB.3p , but not vice versa. In the view of the original problem (|3.1ip , condition (jB.2b ) 
seems to be the natural one. Similarly, in these references the "natural" condition 
(|B.2F ) is replaced by 

P(p,z,t) = 0((/)2 + z2)-i/2) for p2 + ^2 ^ 00, t G M+. (B.4) 

These conditions are in fact equivalent for solutions of (|B.2b ) as becomes clear in the 
subsequent formulation of problem (|B.2p . 

There is an elegant way to eliminate the "coordinate-singularity" at /> = in prob- 
lem ()B.2p . namely by embedding (|B.2p in M^. P is then considered as an axisymmetric 
function in with symmetry axis in xs-direction. Identifying p'^ with X]i=i ^1 
z with X5, X G M^, and introducing Q{x,t) := Q{p,z,t) := P{p, z,t) / the crucial 
observation is [S] 

A,P = {dl --^dp + dl)P = p\dl + -^dp + dl)Q = p^A.Q 

with A5 being the Laplacian in M^. With the further definitions 

bi{x, t) := -Vp{p, z,t) — , i = 1, . . . , 4, b^ix, t) := -Vz{p, z, t). 
P 

c(x, t) := -2 Vp{p, z, t)/p , a{x, t) := r]{p, z, t) 



problem ()B.2p takes in Mp the form 

dtQ-aA5Q = b-VQ + cQ in G5 x M+, (B.5a) 

A5Q = in G5 X M+, (B.5b) 

Q and VQ continuous in x M+, (B.5c) 

Q{x,t) = 0{\x\-^) for |x| ^ 00, t G M+, (B.5d) 

Q{-,0) = Qo , Qo axisym. on G5 x {t = 0}. (B.5e) 



G5 is now an axisymmetric bounded region in M^. An axisymmetric initial field Qo 
implies axisymmetry of Q{-,t) for all t > 0. A condition on the symmetry axis is no 
longer necessary; conditions ()B.2d .e) (as well as ()B.3p ) are automatically satisfied by 
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P := p^Q. As to the behaviour for large x, Q is in G an exterior harmonic function with 
representation (|C.ip : thus, (|B.5l i) imphes |VQ(x, •)! = 0(|x|~'') for |x| — )• oo and hence 
(|B.2f ) as well as ()B.4p . On the other side, in the view of ()C.4p each of the conditions 
dEllf) and dElD implies ^BM ). 

Remark B.l Stredulinsky et al. {16] doubt the correctness of the boundary condi- 
tion ()B.4p in [lOj and cite their own results (theorem 2) about solutions with non- 
vanishing asymptotic value Pooit) at spatial infinity. They mention the possibility of 
limf_>oo Pooit) 7^ 0, even when liuit^Q Poo{t) = 0. In fact, condition (]B.4p is correct as 
demonstrated above, which means Poo = in the axisymmetric problem. The discrep- 
ancy arises because in [16] the authors are especially interested in the two-dimensional 
case where their problem (1) makes physical sense ("dynamo problem with plane sym- 
metry"). In d = 2, in fact. Poo 7^ cannot be avoided in general. In d > 2, however, 
problem (1) is underdetermined and the condition Poo = can be added. Observe in 
this context that in [16] the authors do not claim uniqueness for solutions of problem 
(1), uniqueness is claimed only for weak solutions which are "minimizers" (this is more 
than "harmonic") in G (theorem 1). 



C Exterior harmonic functions 

An exterior harmonic function u is called "harmonic at infinity" iff u{x) — ?■ for 
|x| — 7- oo. In (i > 3 dimensions these functions have the series representation 

oo 

n(x) = J^|x|2-"-'^y„(x/|x|), y„GF„, (C.l) 

n=0 

absolutely and uniformly converging in the exterior of any ball P,. C M'^ with r > 1 
(cf. [SI pi 15]. Hn is the space of all harmonic homogeneous polynomials of degree n 
in M*^ restricted to the unit sphere S"^~^ with dimension P„ := dimP„ = (2n + d — 
2){n + d-3)l[n\{d-2y.]^^ (cf. [HI p98f]. In particular, dimPo = 1 and Yq = const; any 
other Yn has vanishing spherical mean, (1^) = 0, n G N. The total of spaces Hn spans 
L^iS'^-^): L2(5rf-i) = 0^ Hn. So, choosing orthonormal bases {Ynm | 1 < m < P„} 
in Hn, any / G L'^{S'^~^) allows a unique representation 

oo Dn 

f — ^ ^ ^ ^ (^nm Ynm (^'•2) 
n=0 m=l 

with coefficients c„m := if, Ynm) G M. Obviously, ()C.2p is the higher-dimensional 
analogue of the well-known representation 

oo n 

f — ^ ^ ^ ^ Cnm Ynm (^.3) 

n=0 m=—n 

by spherical harmonics {Ynm | ^ G Nq, \m\ < n} in c? = 3. Note that the Ynm are in our 
setting real quantities; by taking suitable linear combinations this is true for spherical 
harmonics as well. 
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Exterior harmonic functions u satisfying the condition Jg^ |Vupdx < oo for some 
R > are harmonic at infinity up to a constant cq (cf. [151 p41])- With (jC.ip this 
imphes the representation 

oo 

n(x) = co + ^|x|2-"-^y„(x/|x|), coeM, Yn€Hn. (C.4) 

n=0 

(]C.4p holds also in the case of exterior harmonic functions with asymptotic conditions 
u{x) = 0(1) for \x\ — )• oo (cf. ^ p64]) or \Vu\ = 0{\x\^^'^) for |x| — )• oo. The latter 
statement follows from the former and the estimate 

\u{x)\ <\u{Rx/\x\)\ + |Vn(r2;/|x|)|dr, \x\>R>0. 



JR 

Let us, finally, consider the exterior boundary-value problem 

Au = in G, (C.5a) 

u = (p on dG, (C.5b) 

u{x) — )• c for l^l — )• oo. (C.5c) 

Here, G C M'^, d > 3 is an exterior region, i.e. G = for some bounded domain G C 

M'^', with C^-boundary dG. For G C(aG) and c G M we call u £ C{G U (9G) n G^{G) 
satisfying ()C.5p a classical solution. A weak formulation is based on the spaces Ti := 
{v G Hl^{G)\\\v\\f^ < oo} and Ho := clos{C7~(G), || • ||^} with \\vf^ := jglVw^dx. 
Describing the boundary and asymptotic conditions by a function w £ Ti, a. weak 
version of (jC.Sp reads: 

/ Vn • Vv dx = for any v G "Hoj (C.6a) 

u-wG-fio. (C.6b) 
For problems ([05]) . holds: 

Theorem C.l (Solution of the exterior Dirichlet problem) The exterior boun- 
dary-value problem (|C.6p with given w £ H has a unique solution u £ H. Moreover, 
u G C°°{G) and Au = in G. If w\qq = </> G G{dG) and w ^ c for \x\ — )• oo, then u is 
a classical solution of ()C.5p . Furthermore, u is the unique minimizer of the functional 
II • 11^ on the set {v + w \ V £ Hq}. 

For a proof we refer to [H p543]. We note only that || • || ^-convergence already im- 
plies II • ||>-2 /^N-convergence, thus Tio C Ti. In fact, the Gagliardo-Nirenberg-Sobolev- 
inequality (cf. [3 p263]) implies the estimate 

II II ^r^ll II - 

for any v G Tio and a constant C depending only on d. So, fixing some X d G we 
obtain: 

\\v\\mK) < C{K)\\v\\Lr>iK) < C(i^)||^||^,(g) < G(K)G||t;||^ . 
In particular, our T-Lq coincides with the corresponding space Bq{G) in [6]. 
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D Poloidal free decay modes 

The poloidal free decay modes are a countable set of explicit solutions of the eigenvalue 
problem (|3.ip if G is a ball Br in M^. In terms of spherical Bessel functions jn and 
spherical harmonics (cf. appendic C) they take the form 



Plr. 




Ynm{x/\x\) m Br 

Jn{ii ) , /GN,nGNo,|m|<n 

(|x|/i?)-"-iy„„(x/|x|) in Br 



with eigenvalues Xinm '■= Mn '■= ^/R)'^'^ is the /-th positive zero of j„. For their 
restrictions pinm ■= PlnmlBa holds 

Theorem D.l The set of functions {pinm ■ Br — ?• M| / G N, n G Nq, |?n.| < n} consti- 
tutes a complete orthonormal system in L^{Br). 

Proof: The orthonormality of the pinm can be checked by explicit calculation using 
the orthonormality of the Ynm and the corresponding relation for the jn (cf. [H p485, 
eq. 11.4.5]. To prove the completeness we show any solution of problem ()3.ip being a 
linear combination of the pinm- Theorem 3.1 yields then the completeness of the pinm- 
So, let u be a solution of (|3.1|) with eigenvalue A > 0. The representation (jC.3p for 
-L^-functions on the unit sphere implies the representation 



oo n 



U{x)='^ Unm{\x\)Ynm{x/\x\) (D.l) 

n=0 171=— n 

for u with coefficients Unm depending only on For |x| > R, (jC.ip yields 

'^nm ( I I ) — (^nm \ X \ 

with Cnm G M, whereas inserting (jP.ip into (|3.1k ) yields the differential equation for 

'^nm on 

1 d / 2 d \ 1 



- :p2-^y' -^'^nmj + :^n{n+l)Unm = >^Unm (D-2) 

with r := After rescaling s := \/Ar, ()D.2p takes the spherical form of Bessel's 
differential equation for Vnm{s) := Unm{s/^/\): 



2 . . 

n(n + l)\vr. 



2^ do 
s^—, + 2s— + s^ 



ds^ ds 

with (in s = 0) regular solutions jn-, n G Nq. Matching inner and outer solutions 
according to (|3.1b ) yields the condition 

(5^Jn(s) + (n + l)j„(s)) [ = (si„-i(s))|^^^^ = (D.3) 



ds" 

fixing the eigenvalue at 



--^/XR 

A = Az„ = {ir'/Rf- (D.4) 
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Note that (|D.3p holds also in the case n = with j_i = cos s/s. 

According to theorem 3.1, the eigenvalue A is of finite multiplicity. Thus, only 
finitely many pairs {l,n) satisfy (jEl), and ([DTT]) is in fact a (finite) linear combination 
oi the pinm- n 

Remark D.2 The condition of vanishing spherical mean eliminates all n = modes; 
thus, {pinm\l G N, n G N, |m| < n} is a complete orthonormal system in {v G 
L^{Br) I (v) = 0}. 

Acknowledgements 

Discussions with and comments from M. Seehafer, C. G. Simader and W. von Wahl 
are gratefully acknowledged. 

References 

[1] M. Abramowitz and I. A. Stegun. Handbook of mathematical functions with for- 
mulas, graphs, and mathematical tables, volume 55. U.S. Government Printing 
Office, Washington, D.C., 1964. 

[2] S. Axler, P. Bourdon, and W. Ramey. Harmonic function theory, volume 137 of 
Graduate Texts in Mathematics. Springer- Verlag, New York, 1992. 

[3] G. Backus. The axisymmetric self-excited fluid dynamo. Astrophys. J., 125:500- 
524, 1957. 

[4] G. Backus. A class of self-sustaining dissipative spherical dynamos. Ann. Physics, 
4:372-447, 1958. 

[5] S. Chandrasekhar. On force-free magnetic fields. Proc. Nat. Acad. Sci. U. S. A., 
42:1-5, 1956. 

[6] R. Dautray and J.-L. Lions. Mathematical analysis and numerical methods for 
science and technology. Vol. 1. Springer- Verlag, Berlin, 1990. 

[7] L. C. Evans. Partial differential equations, volume 19 of Graduate Studies in 
Mathematics. American Mathematical Society, Providence, RI, 1998. 

[8] G. B. Folland. Introduction to partial differential equations. Princeton University 
Press, Princeton, NJ, second edition, 1995. 

[9] D. Gilbarg and N. S. Trudinger. Elliptic partial differential equations of second 
order. Springer- Verlag, Berlin, 1977. Grundlehren der Mathematischen Wis- 
senschaften. Vol. 224. 



27 



[10] D. J. Ivers and R. W. James. Axisymmetric antidynamo theorems in compressible 
nonuniform conducting fluids. Philos. Trans. Roy. Soc. London Ser. A, 312:179- 
218, 1984. 

[11] R. Kaiser. The non-radial velocity theorem revisited. Geophys. Astrophys. Fluid 
Dyn., 101 (3-4): 185-197, 2007. 

[12] D. Lortz and R. Meyer-Spasche. On the decay of symmetric dynamo fields. Math. 
Methods Appl. Sci., 4(l):91-97, 1982. 

[13] H. K. Moffatt. Magnetic Field Generation in Electrically Conducting Fluids. Cam- 
bridge University Press, 1978. 

[14] B. J. Schmitt. The poloidal-toroidal representation of solenoidal fields in spherical 
domains. Analysis, 15(3):257-277, 1995. 

[15] C. G. Simader and H. Sohr. The Dirichlet problem for the Laplacian in bounded 
and unbounded domains, volume 360 of Pitman Research Notes in Mathematics 
Series. Longman, Harlow, 1996. 

[16] E. W. Stredulinsky, R. Meyer-Spasche, and D. Lortz. Asymptotic behavior of 
solutions of certain parabolic problems with space and time dependent coefficients. 
Comm. Pure Appl. Math., 39(2):233-266, 1986. 



28 



